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We have studied the effect of atomic electrons on the nuclear transition from the isomeric 229m Th 
state to the ground 229fl Th state in 229 Th + due to the electronic bridge process. The exact value of 
the nuclear transition frequency is unknown so far; therefore, we have developed a formalism that 
can be used for any nuclear transition frequency. We have calculated positions of several high-lying 
even-parity states which are not presented in experimental atomic spectra databases. We have found 
their energy levels and g factors. 
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The energy splitting of the ground-state doublet of 
the 229 Th nucleus is only several electron volts At 
the same time the exact value of the frequency for the 
transition from the isomeric 229m Th state to the ground 
2299 Th state is unknown. Experiments give values of this 
frequency ujjy ranging from 3.5 ± 1.0 eV Q to 7.6 ± 0.5 
eV 3] . The measurements of the lifetime of the isomeric 
state performed by different experimental groups lead to 
values which differ from each other by many orders of 
magnitude (see, e.g., 0, Q). 

As was noted in Q the nuclear transition from the 
isomeric state to the ground state is of a great interest 
since it makes it possible to build a very precise nuclear 
clock. This transition is very sensitive to hypothetical 
temporal variation of the fundamental constants Q. 

The triply ionized 232 Th was recently laser cooled 0. 
Further, this experimental group plans to investigate the 
nuclear transition between the isomeric and the ground 
state in a cold 229 Th 3+ ion. Another experimental 
group @ plans to use the ion 229 Th + to study the nuclear 

229m Th _229g Th transition _ 

In our previous work ^10| we considered the 229 Th 3+ 
ion and calculated the transition probability of the 229 Th 
nucleus from its lowest energy isomeric state 229m Th to 
the ground state 229 fTh due to the electronic bridge (EB) 
process. In this paper we consider the more complicated 
three- valence ion 229 Th + . In our approach we do not fix 
the value of the nuclear transition frequency. Hence, the 
result obtained here can be applied for any value of uin- 

The paper is organized as follows. In Sec.|TI]we briefly 
discuss the general formalism describing the EB process. 
In Sec. IIIII we describe the method of calculation of the 
properties of Th + . Section IIVI is devoted to the results 
of calculations and Sec. [V] contains concluding remarks. 
Atomic units {% = \e\ = m e = 1 and the speed of light 
c = 137) are used throughout. 



FIG. 1: Feynman diagrams of the EB process. The single 
and double solid lines relate to the electronic and the nuclear 
transitions, correspondingly. The dashed line is the photon 
line. 



II. GENERAL FORMALISM 

A derivation of the equation for the probability of the 
EB process, Teb, is given in detail in [ljj]. For this rea- 
son we will repeat here only the main features of the 
formalism. 

The EB process can be represented by the two Feyn- 
man diagrams in Fig. [TJ In the following we assume that 
the initial i and the final / electronic states are of oppo- 
site parity and fixed. A real photon which is emitted or 
absorbed is the electric dipole photon. The EB process 
can be effectively treated as the electric dipole i — »• f 
transition of the electron accompanied by the nuclear 
transition from its isomeric state to the ground state. 

Because the exact value of the nuclear transition fre- 
quency is unknown we do not fix it in our calculation. 
Using the experimental data we suggest that most prob- 
ably the real value of lun is between 2 and 8 eV. The 
general expression for Teb we used for calculation of the 
EB process for 229 Th 3+ can be simplified for 229 Th + . 
This is due to the spectrum of Th + being much denser 
than the spectrum of Th 3+ . As a result, for any nuclear 
transition frequency wjv lying between 2 and 8 eV we 
can find an atomic transition from the initial state i to 
the definite intermediate state n whose frequency will be 
very close to a; at. Assuming the resonance character of 
the EB process we arrive at the following expression for 
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where .Mi is the magnetic dipole nuclear moment and 
\I g ) and \I m ) are the ground nuclear state and the iso- 
meric nuclear state (I g = 5/2 + , [633] Nilsson state and 
Im = 3/2 + , [631] Nilsson state); Ji is the electron total 
angular momentum of the initial state, w is the real pho- 
ton frequency determined from the law of conservation 
of energy as u> = £j — Sf + uin (where Ek is the atomic 
energy), and G\ can be approximated by 
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Here 71 is the electronic magnetic-dipole hyperfme cou- 
pling operator. The total hyperfine coupling Hamiltonian 
Hhfi may be represented as 
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The operator D is the electric dipole moment operator, 
ujin = Si — e n , and 7fe encapsulates all other electronic 
quantum numbers. The explicit expressions for the ma- 
trix elements of the operators T\ and D are given in our 
paper (To| . If we introduce the quantity 



R n = K-ffJfWDW^j^^JnWTiHJJl 2 

then Eq. ([2]) can be rewritten as 
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In Eqs. ©-© the electronic state \j n J n ) is assumed to 
be fixed. This state should be chosen to meet two con- 
ditions: 1) —u>in ~ ujn and 2) if the first condition is 
fulfilled for two atomic states we should take the state 
for which R n is larger. The second condition is impor- 
tant because in certain cases the coefficients R n for two 
neighboring energy levels differ by several orders of mag- 
nitude. 

In Ref. [13] we used the dimensionless quantity /3 M1 
defined as the ratio of the probability of the EB process, 
Teb, to the probability of the Ml radiative nuclear m — > 
g transition, T^: 
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It is reasonable to choose the ground state (6cP 7s) 
J = 3/2 as the initial state i and consider the lowest lying 
odd-parity state (5/ 7s 2 ) J = 5/2 as the final state /. 
Thus, the intermediate atomic states contributing to Gi 
are even-parity states and our purpose is to calculate the 
coefficients R n for all even-parity states whose transition 
frequencies to the ground state are between 2 and 8 eV. 
Then, using Eqs. ([5]) and ([B]) we can find the quantities 
Gi and ft mi, correspondingly, for any w^r lying between 
2 and 8 eV. 



III. METHOD OF CALCULATION 

We consider Th + as the atom with three valence elec- 
trons above closed-shell core [Is 2 , . . . ,6p 6 ]. We employ 
the approach combining the configuration-interaction 
(CI) method in the valence space with many-body per- 
turbation theory (MBPT) for core polarization effects. 
In the following we refer to this combined approach as 
the CI+MBPT method [HI- 

At the first stage we have solved Dirac-Hartree- 
Fock (DHF) equations [l2| in the V N ~ 3 approximation. 
This means that the DHF equations were solved self- 
consistently for the core electrons. After that we de- 
termined the 5/, 6d, Ip, 7s, and 8s orbitals from the 
frozen-core DHF equations. The virtual orbitals were 
determined with the help of a recurrent procedure [l3j . 
The one-electron basis set included ls-18s, 2p-17p, 3d- 
16d, and 4/-15/ orbitals on the CI stage. 

The configuration spaces for even-parity and odd- 
parity states were formed as follows. The main config- 
uration of the ground state is 6d 2 7s. We formed the 
configuration space for the even-parity states by allow- 
ing single, double, and triple excitations from the 6d 2 7s 
configuration to the 7s-13s, 7p-12p, 6d-lld, and 5/-10/ 
shells. The main configuration of the lowest lying odd- 
parity state is 5/7s 2 . The configuration space for the 
odd-parity levels was formed by single, double, and triple 
excitations from the 5/7s 2 configuration to the 7s-13s, 
7p-12p, 6c?— lid, and 5/-10/ shells. Inclusion of all possi- 
ble (up to triple) excitations is important, especially for 
high-lying states. It allows us to take into account most 
completely the configuration interaction for all consid- 
ered states. 

In the CI+MBPT method, the energies and the wave 
functions are determined from the eigenvalue equation in 
the model space of the valence electrons, 

tf cff (£ p )|<& P H£ p |<I> p ), (7) 

where the effective Hamiltonian is defined as 



H< B (E) = H FC + E(E). 



(8) 



Here H^c is the relativistic three-electron Hamiltonian 
in the frozen-core approximation and £(£?) is the energy- 
dependent core-polarization correction. 

Together with the effective Hamiltonian H r we intro- 
duce the effective electric-dipole operator D c g and the 
operator (7i) e ff acting in the model space of valence elec- 
trons. These operators were obtained within the rel- 
ativistic random-phase approximation (RPA) 0, [l5| . 
which describes a shielding of the externally applied elec- 
tric field by the core electrons. The RPA sequence of 
diagrams was summed to all orders of the perturbation 
theory. 

To solve the RPA equations and to calculate diagrams 
for the effective Hamiltonian and the effective operators 
D and 71 we used a different basis set. The core orbitals 
in this basis set are the same as before, but the number of 
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virtual orbitals is much larger. On the whole, it consisted 
of ls-22s, 2p-22p, 3d-22d, 4/-22/, and 5g-16g orbitals. 

IV. RESULTS AND DISCUSSION 

We start the discussion of the results with the following 
remark: The spectrum of Th + is very complicated. As is 
seen from the experimental data [la ], on the one hand, 
the states belonging to different configurations strongly 
interact with each other and LS coupling is not valid 
(even approximately) for this ion. On the other hand, 
it is not a chaotic system. Respectively, the methods 
of statistical physics are not applicable. Such an "inter- 
mediate" type of coupling makes the calculations of the 
properties of Th + rather difficult. 

As we have already mentioned in Sec. HU we consider 

the following transition: Qd 2 7s{J = 3/2) — i> n -^4 
5f7s 2 (J = 5/2). According to Eqs. © - © only inter- 
mediate states n with J n — 3/2 and J n = 5/2 contribute 
to the probability of the EB process for this transitions. 

In Tables |T] and [TT| we presented the calculated values 
of the energy levels with J n = 3/2 and J n — 5/2 and 
also g factors and the coefficients R n obtained with use 
of Eq. © for the most interesting frequency range from 
2 to 8 eV. In Table UJ we present the results for the atomic 
frequencies from 2 to 5 eV and in Table [TT| (which is a 
continuation of Table U} the data are listed for the fre- 
quencies from 5 to 8 eV. The results for the energy levels 
and g factors were obtained in the CI+MBPT approxi- 
mation. The values of the coefficients R n were found in 
the frame of the CI+MBPT+RPA approach. 

As is seen from the tables basically the agreement be- 
tween the experimental and the calculated energy levels 
is satisfactory. For the majority of the levels presented 
in Tables HI and [TTI the agreement is at the level of several 
percent. The largest difference between the experimental 
and the theoretical values is for the states belonging to 
the 6d 3 configuration, where it reaches 15%. At the same 
time the g factors for these states were reproduced rather 
well. This means that the configuration interaction was 
taken into account correctly. 

The energy levels with total angular momenta J = 3/2 
and J = 5/2 lying higher than 56391 cm -1 are not iden- 
tified experimentally. In our work we have determined 
several new high- lying energy levels with J = 3/2 and 
5/2. In the first rows of Tables HI and [TTI we indicate the 
configurations that give the largest contributions to these 
states according to our calculation. 

As we have already mentioned the configuration mix- 
ture is strong for all states starting from the ground state. 
Sometimes we were unable to reproduce correctly the 
configuration interaction. In such cases the theoretical g 
factors differ from the experimental g factors and, respec- 
tively, the accuracy of calculation of R n for such states 
is poorer. 

As follows from Tables UJ and [TTJ the coefficients R n 
change from 10 -7 to 10 -1 . This is not surprisingly if 



TABLE I: The low-lying energy levels in the range from 18119 
to 40644 cm" 1 (from 2 to 5 eV) in the CI+MBPT approxi- 
mation, g factors, and the coefficients R n (in a.u.). A is the 
difference between the energies of the ground state and the ex- 
cited state. The notation y[x] means y x 10 s . The theoretical 
values are compared with the experimental data. 
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we note that the initial state 6d 2 7s and the final state 
5/7s 2 differ from each other by two electrons while T\ 
and D are the one-electron operators. For this reason 
the i — > n — > f transition occurs only by the configura- 
tion interaction. In a case when the intermediate state n 
is characterized by configurations that open two strong 
one-electron 6d 2 7s — > n and n — > 5/7s 2 transitions, R n 
turn out to be large. Due to the complexity of the energy 
level spectrum of Th + the accuracy of the calculation of 
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TABLE II: The low-lying energy levels in the range from 40924 to 64000 cm" 1 (from 5 to 8 eV) in the CI+MBPT approximation, 
g factors, and the coefficients Rn (in a.u.). A is the difference between the energies of the ground state and the excited state. 
The notation y[x] means y x 10^. 
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the coefficients R n is not high. We would consider these 
values as an order-of-magnitude estimate. 

To illustrate how the developed formalism works we 
consider two possible values of the nuclear frequency, 
ujn = 3.5 eV HI and ujn = 5.5 eV [Vq , as reported by 
two experimental groups in the mentioned papers. In 
Table ITO1 we present the values of the relevant quantities. 

For wjv = 3.5 eV « 28231 cm -1 the resonance con- 
tribution to Teb comes from the atomic state J = 5/2 
at 28824 cm -1 belonging to the configuration 5f7s7p. 
We chose this state because the transition frequency lo ics 
from this state to the initial state i (the ground state) is 
close to u>n and the coefficient R n is largest. Knowing 
from Table Q] the coefficient R n for this state and using 
Eqs. ([5]) and ((HJ we can easily find the quantities G\ and 

/3 M i for the transition 6d 2 7s (J = 3/2) -A- 5f7s7p (J = 

5/2) 5/7s 2 (J = 5/2). In a similar way Gi and ft mi 
can be obtained for ujn — 5.5 eV. 

Comparing the coefficients (3mi obtained for lun — 
3.5 eV and ujn = 5.5 eV we see that they are of the 
order of 10 2 -10 3 . We note that in the case of lun = 
5.5 eV the difference (w rcs — uin) is only 26 cm -1 while 
R n = 0.001 is rather small. For ojn — 3.5 eV the 
difference (w res ~ w/v) ~ 600 cm -1 but the coefficient 
R n = 0.2 is two orders of magnitude larger than that for 
ujn = 3.5 eV. The latter occurs because the resonance 
energy level whose frequency is close to ujn = 3.5 eV 
belongs to the configuration 5/7s7p. Hence, there is a 

strong 5/7s7p (J = 5/2) 5/7s 2 (J = 5/2) transition. 
Due to an admixture of the configuration 6c? 2 7s to the 
configuration 5f7s7p the amplitude of the 6d 2 7s (J — 

3/2) — 5f7s7p(J = 5/2) transition is not small. As a 
result, the coefficient R n is large. 

The case of ujn = 7.6 eV [3[ requires special atten- 
tion. The problem is that the atomic energy levels are 



not identified experimentally in the region of 7.5 eV and, 
consequently, we cannot compare the theoretical energy 
levels with the experimental energy levels. As we previ- 
ously mentioned, the theoretical accuracy is at the level 
of several percent. Thus at present we are unable to re- 
liably predict the position of the resonance energy level 
and, consequently, the coefficient (3mi- For this reason 
experimental investigations and identification of the en- 
ergy levels in the frequency region ~ 7.5 eV would be 
very useful. Once these tasks are completed the coeffi- 
cient Pmi can be easily determined. 



V. CONCLUSION 

To conclude, we have found several high- lying even- 
parity states with total angular momenta J = 3/2 and 
J = 5/2 that are not identified in the atomic spectra 
database [l6| ■ We have determined the energy levels and 
the g factors of these states. 

We have calculated the coefficients R n determined by 
Eq. (0J for the even-parity states lying between 2 and 
8 eV. When the nuclear transition frequency ujn is ex- 
actly known and the atomic energy levels are experimen- 
tally identified, we can find, using R n , the coefficients 
G\, Ami, and the probability of the EB process. 
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TABLE III: The nuclear transition frequency ujn (given in eV and in cm -1 ) along with the configuration, the total angular 
momentum J, and the transition frequency with respect to the ground state (u> r es) for the resonance state mainly contributing 
to Gi, listed with the coefficients R n (in a.u.), Gi (in a.u.), and flui- 









Resonance state 










eV 


cm -1 


Conf. 


J 


w res (cm _1 ) 


Rn 


Gi 


Pmi 


3.5 


28231 


5/7s7p 


5/2 


28824 


0.2 


4570 


225 


5.5 


44363 


5/ 2 6d + 5/6d7p 


5/2 


44389 


0.001 


11880 


720 



